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Vibrational vs. electronic coherences in 2D spectrum of molecular systems 
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Two-dimensional spectroscopy has recently revealed the oscillatory behavior of the excitation dynamics of 
molecular systems. However, in the majority of cases there is considerable debate over what is actually being 
observed: excitonic or vibrational wavepacket motion or evidence of quantum transport. In this letter we 
present a method for distinguishing between vibrational and excitonic wavepacket motion, based on the phase 
and amplitude relationships of oscillations of distinct peaks as revealed through a fundamental analysis of the 
two-dimensional spectra of two representative systems. 



Two-dimensional photon-echo (2DPE) spectroscopy is 
a powerful tool capable of resolving quantum correlations 
on the femtosecond timescalei"— . They appear as beats 
of specific peaks in the 2DPE spectrum for a number of 
molecular systems 3,4 . However, the underlying processes 
are often ambiguous. At first, the beats were attributed 
to the wave-like quantum transport with quantum coher- 
ences being responsible for an ultra-efficient excitation 
transfer 3 6 . The same process was associated with the 
opposite phase beats in the spectral regions which are 
symmetric with respect to the diagonal linei 

In molecules and their aggregates, electronic transi- 
tions are coupled to various intra- and intermolecular vi- 
brational modes. Vibrational energies of these are of the 
order of 100 - 3000 cm -1 , while the magnitudes of the res- 
onant couplings, J, in excitonic aggregates (e.g. in photo- 
synthetic pigment-protein complexes or in J-aggregates) 
are in the same range. Thus, vibronic and excitonic sys- 
tems show considerable spectroscopic similarities, and 
presence of electronic and/or vibrational beats in the 
2DPE spectrum is expected. Indeed, similar spectral 
beats originating entirely from a high-energy vibrational 
wavepacket motion have been observed^. The possibil- 
ity of distinguishing the electronic and vibrational origin 
of the beats from a 2DPE spectrum has been empha- 
sized in a recent letter—. However, the reported conclu- 
sions have not been supported by theoretical arguments, 
and thus are questionable. Therefore, the highly rele- 
vant question of how vibrations interfere with electronic 
coherences in 2DPE spectrum is still an open one. A the- 
oretical study of the origin of spectral beats, their phase 
relationships in the rephasing and non-rephasing compo- 
nents of the 2DPE spectrum is presented in this article. 

We address this problem by considering two generic 
model systems which exhibit distinct internal coherent 
dynamics. The simplest model of an isolated molecular 
electronic excitation is the vibronic system represented 
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FIG. 1. Energy level structure of the displaced oscillator (a) 
and electronic dimer (b) and corresponding linear absorption 
spectra. 



by two electronic states, \g) and |e), which are coupled to 
a one-dimensional nuclear coordinate q. We denote the 
model by a displaced oscillator (DO) system (Fig. [TJl) . 
Taking h =1, the vibronic potential energy surface of 
the \e) state is shifted up by electronic transition energy 
uo eg and its minimum is shifted by d with respect to the 
ground state \g); d is the dimensionless displacement. 
This setup results in two vibrational ladders of quantum 
sub-states \g m ) and |e n ), m,n = O...00, characterized 
by the Huang-Rhys (HR) factor HR = d 2 /2±m. 

The other model system, which shows similar spec- 
troscopic properties but has completely different coher- 
ent internal dynamics without vibrations, is an excitonic 
dimer (ED). It consists of two two- level chromophores 
(sites) with identical transition energies e. The two sites 
are coupled by the inter-site resonance coupling J. As 
a result, the ED has one ground state \g), two single- 
exciton states fa) and fa) with energies e eu e 2 = e ± J, 
respectively, and a single double-exciton state |/) with 
energy Sf = 2e — A, where A is the bi-exciton binding 
energy (Fig. [TJd)^. 

The absorption spectrum of both systems is as fol- 
lows. The absorption of the DO is determined by transi- 
tions from the \g m ) vibrational ladder into \e n ) scaled 
by the Franck- Condon (FC) vibrational wavefunction 
overlap s 1 ^ 13 . Choosing HR = 0.3 and k^T « \ujq and 
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assuming Lorentzian lineshapes with linewidth 7, we get 
the vibrational progression in the absorption spectrum 
(dashed line in Fig. [I]). Here uo is the vibrational energy. 
The most intensive peaks at uj eg and uo eg +ujo correspond 
to 0-0 and 0-1 vibronic transitions. Qualitatively similar 
peak structure is featured in the absorption of ED, where 
the spectrum shows two optical transitions \g) — >• \e\) 
and \g) — >• |e2), assuming both are allowed. Choosing 
J = ujq/2 and the angle (p between the chromophore tran- 
sition dipoles equal to 7r/6, and using adequate linewidth 
parameters, we get absorption peaks (solid line in Fig. [1]) 
that exactly match the strongest peaks of the DO. As 
expected one cannot distinguish between these two inter- 
nally different systems from the absorption spectra alone. 

The 2DPE spectrum carries more information than ab- 
sorption. However, it consists of many contributions and 
unambiguous distinction between the ED and DO sys- 
tems becomes difficult. In order to unravel the 2DPE 
spectra we thus need to construct the entire 2D signal 
from the first principles for both systems and recover the 
source of oscillations in the 2DPE spectrum. 

In the conventional scheme of the 2DPE measurement, 
two primary excitation pulses with wavevectors k\ and 
&2 followed by the probe pulse ks are used; kj are pulse 
wavevectors. The signal is detected at the ks = —fci + 
&2 + &3 phase-matching direction. The order of fci and &2 
defines the rephasing configuration (fci) when ki comes 
first and the non-rephasing configuration (fen) when &2 
comes first. 

Semiclassical perturbation theory with respect to the 
incoming fields reveals the system-field interaction and 
evolution sequences, often denoted by the Liouville space 
pathways. Three types of distinct interaction config- 
urations are denoted by the Excited State Emission 
(ESE), Ground State Bleaching (GSB) and Excited 
State Absorption (ESA) contributions^^. If we ne- 
glect environment-induced relaxation, the signals are 
given as sums of resonant contributions, 5(^3,^2,^1) = 
En ^n(^3, £2, u\) of the type 
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FIG. 2. (a) Scheme of contributions to 2DPE spectrum of 
the fei and fen signals for the reduced DO and ED (A = 
0) systems. The ESE contribution is indicated by squares, 
GSB - diamonds, ESA - circles. Solid symbols denote non- 
oscillating contributions in £2, open - oscillatory in the form 
of ±cos(£2£2), where £2 = cjo f° r DO and £2 = 2 J for ED. 
(b) Phase <j> of the contribution (Eq. [3j) and peak profile S n 
as a function of the shift from the peak center (si = S3 = 0) 
using relative coordinates. The diagonal lines of the fei and 
fen contributions to the 2D spectra are shown by dashed lines. 
Peaks are labeled in plots as '1-1', '1-2', etc. 



S n (u3,t 2 , Ul ) = A (n) J J dt 1 dt 3 e +iW3t3+iwitl 

x[±G 3 (t3)G2(t2)Gl(tl)]( n)) (1) 

where the subscript n denotes different terms of the sum- 
mation. A(ji) is a complex prefactor, given by the tran- 
sition dipoles and excitation fields, the propagator of the 
density matrix G for the jtfi (j = 1,2,3) time delay is of 
the one-sided exponential function type 



Gjitj) = 0(tj) exp(-ie jtj) 



(2) 



£j coincides 



(#(£) is the Heaviside step- function). Here 
with the energy gap cj a & between the left and right states 
of the system density matrix relevant to the time interval 
tj. ESE and GSB carry sign while ESA has '-' overall 
sign. 

The Fourier transforms in Eq. (pQ) map the con- 
tributions to the frequency- frequency plot (^,£3) — » 



(uji^ujs) ~ (Tkil^s) (the upper sign is for fci, the lower 
- for fen). Diagonal peaks at uo\ = ^ujs are usually 
distinguished, while the anti-diagonal line is defined as 
=fui + us = Const. The whole 2DPE signal becomes 
a function of £2: either oscillatory for density matrix 
coherences \a)(b\ with characteristic oscillation energy 
£2 = ^ab 7^ 0> or static for populations \a)(a\ (£2 = 0). 

To reveal oscillatory contributions in the DO and ED 
systems we have grouped all contributions into either 
oscillatory or static as shown in Fig. [2^. If we con- 
sider only the two main vibrational sub-states in DO, 
the 2DPE spectrum will have only ESE and GSB contri- 
butions, while ED additionally has ESA. As a function 
of £2, the DO system has 8 oscillatory and 8 static con- 
figurations, which organize into six peaks, while the ED 
system has only 4 oscillatory and 8 static contributions 
which give four peaks. The net result is that the diagonal 
peaks in the fei and cross-peaks in the fen signals are non- 
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oscillatory in the ED, while all peaks except the upper 
diagonal peak in k\ are oscillatory in the DO. We thus 
find significant differences in oscillatory peaks between 
ED and DO systems. 

An important additional parameter to consider is a 
phase of oscillation. Eq. ([1]) can be analytically inte- 
grated for Gj(tj) oc exp(— iSjtj — jjtj). For a single con- 
tribution S n giving rise to a peak at (cji,^) = (=f|£i|> £3) 
we shift the origin of (cji,^) plot to the peak cen- 
ter by introducing the displacements (uji + e\ = — si, 
^3 — £3 = S3 for the rephasing pathways, while uji — S\ = 
si, 0J3 — £3 = S3 for the nonrephasing). For 7 « 71 « 73 
we get the peak profile 

S n (s3,t 2 , si) = A n L(si,s 3 )e" 72£2 cos(|£ 2 |*2 + 0(si,s 3 )), 

(3) 

where the lineshape and phase for the fei (upper sign) 
and k\\ (lower sign) signals are 



L(si,s 3 ) 



V[7 2 ±^3] 2 +7 2 feT^i) 2 



(*?+7 2 )(*|+7 2 ) 



^( s i> s 3) = sg n { £ 2) arctan 



/ 7(s3 T 



V(TSiS3 - 7 2 ) 



(4) 
(5) 



The phase <p and the full profile for A n = 1 and £2=0 
are shown in Fig. [2k. The rephasing and non-rephasing 
configurations are obtained by flipping the direction of 
the s\ axis. At the center of the peak (s\ = S3 = 0), we 
have and = 0, leading to S n oc cos(|£2|^2)- However, 
for (si ^ 0,s 3 ^ 0) we find S n oc cos(|^ 2 | ^2 + 0(si,s 3 )) 
with 0(si, S3) 7^ 0. Thus, £/ie displacement from the peak 
center determines the phase of the spectral oscillations. 
Note that the sign of the phase <p is opposite for the peaks 
above (£2 < 0) and below (£2 > 0) the diagonal line, and 
this applies for all contributions. 

The whole 2DPE spectrum is a sum of all relevant 
contributions. Assuming that all dephasings are similar, 
different contributions to the same peak will have the 
same spectral shape and they may be summed. We can 
then simplify the 2DPE plot by writing the signal as a 
sum of peaks which have static (from populations) 
and oscillatory (from coherences) parts: 



S(L>3,t2,U>l) 



-72*2 



/ Lij((Ji,U)3) 



A%j • cos 



Here 



(\Uij\t2 + (j>ij(ui,U)3))] . (6) 

is the characteristic oscillatory frequency of 
a peak (ij), A^fe) and A^fa) are the real parts of 
orientationally-averaged prefactors of population and co- 
herence (electronic or vibronic) contributions, respec- 
tively. The spectral lineshape is given by Lij((Ji, 0J3). 
Here we clearly identify the oscillatory amplitude and its 
phase for a specific peak. 

To apply this expression to our systems, we assume a 
typical situation where the spectrum of the laser pulses is 
tuned to the center of the absorption spectrum and the 
limited bandwidth selects the two strongest absorption 
peaks. In the 2DPE spectra two diagonal and two off- 
diagonal peaks for ED and DO are observed. Indices i 
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FIG. 3. The amplitudes of oscillatory peaks of 2D spectra of 
the DO model for k\ and k\\ signal. Note that the negative 
amplitude denotes a phase shift of 7rof the oscillation. 



and j in Eq. (|6j) run over the positions of the peaks and 
thus can be (1,1), (1,2), (2,1), and (2,2). For clarity we 
study the spectral dynamics with £2 at the short delays, 
t2 <C 72 -1 , and use notations A, L for the k\ signal and 
A, L for the ku signal. 

The transition dipole properties of the ED results in 
the picture where all static amplitudes of the ED are posi- 
tive and A^ = A^ , A22 = A^ 2 •> ^-12 = ^21 = ^21 = A^ 2 • 
The oscillatory amplitudes are equal: A\ 2 = ^21 = 
A\ 1 = A\ 2 . Such relationships are obtained by consider- 
ing the all-parallel organization of polarization of incom- 
ing electric fields and neglecting the bi-exciton binding 
energy. The spectral beats with £2 can thus only have 
the same phases in the k\ or fen spectrum, when mea- 
sured at peak centers. Additionally the oscillatory ESE 
and ESA parts in ED cancel each other if A = and their 
broadenings are equal. As these relationships do not de- 
pend on coupling J and transition dipole orientations, 
all ED systems should behave similarly. By studying the 
whole parameter space, it can also be shown that these 
relations hold for a hetero-dimer. 

The amplitude-relationships, however, are different for 



the DO system. The amplitudes A\^ of the oscillatory 
peaks are plotted in Fig. [3] as a function of the HR factor, 
where now we include all vibrational levels in the \g m ) 
and \e n ) ladders. For fei, the amplitudes A\ x and A 22 
maintain the opposite sign when HR < 2 and are both 
positive when 2 < HR < 3 (note that A 22 = when 
only two vibrational levels are considered in Fig. (2Jd). 
The oscillation amplitudes A\ x and A 22 change sign at 
HR — 1. Amplitudes A\ 2 and A 21 are always positive. 
Spectrum oscillations with £2 f° r both diagonal peaks in 
the fen signal will be in-phase for the whole range of the 
HR factor. The same pattern holds for the 1-2 and 2-1 
cross-peaks, which will oscillate in-phase, but will be of 
opposite phase compared to the diagonal peaks in the 
region of HR < 1. Note that the sign of amplitudes 
changes with the HR factor, since the overlap integral 
between vibrational wavefunctions can be both positive 
and negative. The amplitudes of static contributions are 
positive in the whole range of parameters and are iden- 
tical for both fei and fen signals. 
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FIG. 4. 2DPE spectra and peak values of the DO and ED as 
the function of population time £2, of the k\ and fen signals. 
Spectra are normalized to the maxima of the total spectra of 
the DO and ED. 



We thus find very different behaviour of oscillatory 
peaks of DO and ED systems. The above analysis ap- 
plies for the central positions of the peaks, which may 
be difficult to determine if the broadening is large. Note 
that the phase <\> varies from — tt/2 to +7r/2 (Eq. ((3]) and 
Fig. [2h) when probing in the vicinity of the peak. How- 
ever, = along the diagonal line for fei and along the 
anti-diagonal line for fen. These lines can thus be used 
as guidelines for reading phase relations of distinct peaks 
in the 2DPE spectrum. For instance, the two diagonal 
peaks can be calibrated by reading their amplitudes at 
the diagonal line, or the two opposite cross-peaks can be 
compared by drawing anti-diagonal lines. 

The 2DPE spectra for both DO and ED systems cal- 
culated by including phenomenological relaxation and 
Gaussian laser pulse shapes 14 i 15 are plotted in Fig. |4j 
The structure and the £2 evolution of the spectra illus- 
trate the dynamics discussed above and clearly shows the 
distinctive spectral properties of the vibronic vs. elec- 
tronic system: (i) diagonal peaks in the fei signal are 
oscillating in DO, but only exponentially decaying in ED 
(the oscillatory traces come from the overlapping tails 
of off-diagonal peaks), (ii) the relative amplitude of os- 
cillations is much stronger in DO as compared to ED, 
where the ESA and ESE cancellation suppresses the os- 
cillations, (iii) opposite oscillation phases are observed in 
DO, while all peaks oscillate in-phase in ED. 

The up-to-date experiments are capable of creating 
broad-band pulses^. Thus, the overtones in DO can be 
excited and beats of nuo frequencies (n is integer) ob- 



served. These may become important in the case of large 
HR factors. Such frequencies are absent in the ED sys- 
tem, since only one oscillatory frequency, equal to 2 J is 
available. 

The analysis presented in this article provides a clear 
physical picture of electronic and vibronic coherence 
beatings in 2DPE spectra. We are able to discrimi- 
nate weakly damped electronic and vibronic coherent 
wavepackets in molecular systems based on fundamental 
theoretical considerations. Dynamics of diagonal peaks 
and cross-peaks as well as relative phase between them 
in the rephasing signal can now be classified for vibra- 
tional and excitonic systems as follows, (i) Static diago- 
nal peaks and oscillatory off-diagonal peaks signify pure 
electronic coherences, not involved in energy transport, 
(ii) Oscillatory diagonal peaks in accord with off-diagonal 
peaks (0 or tt phase relationships) signify vibronic origins. 
The oscillation phase is for electronic coherences and 
or 7r for vibronic coherences. These outcomes hold if the 
signal is probed at the very centers of the spectral reso- 
nances. The observed phase of the beatings varies as the 
signal is recorded away from the center of an oscillating 
peak. 

Our results might be useful in analysis of recently ob- 
served beatings in molecular systems. For instance phase 
relations of the beatings detected at separate points in 
the vicinity of the same cross-peak of the photosynthetic 
LH2 complex 16 might be the result of measurement away 
from the peak center (see Fig. [2}d). The issue of probing 
away from the peak centers also applies to the opposite- 
phase beatings reported by Collini et al. 7 . Our analysis 
thus shows that the detailed phase relationships in the 
two dimensional spectra may be of critical importance. 
By helping to identify spectral beats in photosynthetic 
aggregates, the presented analysis should facilitate an- 
swering the question of importance of electronic coher- 
ences in excitonic energy transfer, its efficiency and ro- 
bustness. 
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